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Abstract. Cells moving on a 2dimensional substrate generate motion by 
polymerizing actin filament networks inside a flat membrane protrusion. New filaments 
are generated by branching off existing ones, giving rise to branched network structures. 
We investigate the force-extension relation of branched filaments, grafted on an elastic 
structure at one end and pushing with the free ends against the leading edge cell 
membrane. Single filaments are modeled as worm-like chains, whose thermal bending 
fluctuations are restricted by the leading edge cell membrane, resulting in an effective 
force. Branching can increase the stiffness considerably; however the effect depends on 
branch point position and filament orientation, being most pronounced for intermediate 
tilt angles and intermediate branch point positions. We describe filament networks 
without cross-linkers to focus on the effect of branching. We use randomly positioned 
branch points, as generated in the process of treadmilling, and orientation distributions 
as measured in lamellipodia. These networks reproduce both the weak and strong 
force response of lamellipodia as measured in force-velocity experiments. We compare 
properties of branched and unbranched networks. The ratio of the network average of 
the force per branched filament to the average force per unbranched filament depends 
on the orientation distribution of the filaments. The ratio exhibits compression 
dependence and may go up to about 4.5 in networks with a narrow orientation 
distribution. With orientation distributions measured in lamellipodia, it is about 2 
and essentially independent from network compression, graft elasticity and filament 
persistence length. 
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The crawling of many different cell types is essential for life. Undifferentiated cells move 
towards the site, where they form a tissne or organ in the developing embryo. Skin 
cells start crawling when they have to close a wonnd [1]. Dnring metastasis, cancer cells 
dissociate from the primary tnmor, crawl towards blood vessels and spread all over the 
body [2, 3]. Branched actin hlaments carry forces dnring cell motion, and conseqnently 
nnderstanding their elastic properties is central to nnderstanding the mechanics of cell 
motility. 

In vitro, cells are plated on a two dimensional snbstrate to observe their dynamics. 
They form a flat membrane protrnsion in the direction of motion, the lamellipodinm, 
which is only abont 100-200 nm thick bnt several pm deep and wide [4]. A dense 
network of actin hlaments (F-actin) inside the lamellipodinm pushes the leading edge 
membrane forward [5]. Treadmilling of the hlaments drives motion [6]: Filament barbed 
(or plus) ends polymerize at the leading edge of the lamellipodinm and the pointed (or 
minus) ends depolymerize at the rear. 

Usually cells move in response to an external signal. A variety of signals stimulate 
the activation of nucleation promoting factors (NPFs) (like WASp or WAVE) located in 
the leading edge membrane of the lamellipodinm. They activate the actin related protein 
complex Arp2/3. It binds to an existing hlament very close to or at its barbed end at 
the lamellipodium’s leading edge. That initiates the growth of a new hlament branch 
out of the Arp2/3 complex. Many of these branched structures consisting of mother 
hlament and branch form the F-actin network in the lamellipodinm. The branched 
structure itself is dynamic. The branch point with the Arp2/3 complex moves rearward 
due to treadmilling in the same degree as mother hlament and branch grow. Since 
Arp2/3 binding to the individual hlaments is not synchronous we hnd at any time many 
diherent positions of branch points in the lamellipodinm F-actin network. 

The elastic properties of the F-actin network crucially depend on the density of links 
between hlaments [7, 8, 9, 10, 11, 12, 13]. Molecular links arise in two ways: Cross¬ 
links connecting two hlaments at some point along their contour length are formed by 
cross-linker molecules like hlamin or a-actinin, and branching attaches the minus end 
of a hlament laterally to a mother hlament. Intuition suggests that branching alone 
could stihen the network to some degree, since branching is a geometrical constraint 
on the conhguration of two hlaments. That intuition has never been quantihed before, 
but is supported by our results presented below. On the other hand, the network 
region close to the leading edge was found to be as soft as weakly cross-linked actin 
networks [9, 10, 11, 12, 13], and experiments in actin solutions suggest that branching 
contributes very little to the elastic modulus of F-actin networks [14]. Here, we would 
like to present a hrst step in quantifying the contribution of branching to the elastic 
and semi-hexible properties of the lamellipodial F-actin network. How much stiher than 
single hlaments are branched hlaments? How are their properties rehected in network 
behavior? We will answer these questions by investigating a single branched hlament and 
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networks of branched structures in an approximation neglecting interactions between 
them in order to focus on branching effects. This neglect of interactions implies that we 
consider only elastic properties on short time scale and not the visco-elastic properties 
arising from cross-linking. 

The mother hlament is grafted at one end and has a free tip at the other one 
in our model system (see Fig. 1). The graft is provided by a highly cross-linked part 
of the F-actin network. This idea is based on the increasing degree of cross-linking 
and hlament bundling towards the rear of the lamellipodium, which has been observed 
in many different experiments and simulations [5, 15, 16, 17, 18, 19, 20, 21, 22, 23] 
(see [24] for a detailed discussion). The graft moves in the direction of cell motion due 
to cross-linker binding and bundling, and thus contour length ’hows’ into the graft. In 
the steadily moving cell, the balance between polymerization and cross-linking creates 
a stationary distance between graft and leading edge membrane. At any time, the 
branched structure is in a conhguration similar to Fig. 1, but the branch point moves 
with treadmilling towards the graft point. Both mother hlament and branch polymerize 
such that their barbed ends stay at the leading edge membrane. Consequently, we will 
vary the position of the branch point on the mother hlament from 0 to its full length 
L when we investigate the elastic properties. We model the hlaments in the weakly 
bending regime, i.e. bending does not ahect the end-to-end distance. We allow for 
an elastic graft of stihness Kg and model the membrane by a constraint, enforcing the 
hlaments to be entirely on the left side of the membrane. 

After the introduction of the model, we will discuss results for a single branched 
hlament. The F-actin network in the lamellipodium comprises hlaments with many 
diherent tilt angles and branch point positions. Hence, we consider it as an ensemble of 
branched structures with varying branch point and tilt angle and calculate its properties 
as the average over this ensemble. The network average of the force will be compared 
to experimental results from force-velocity measurements and the force of unbranched 
networks. 

2. Results 

2.1. The model: semiflexible branched filaments in the weakly bending limit 

We model the interior of the lamellipodium as a two-dimensional space, since it is 
approximately hat as described above. Furthermore, the results are easily generalised 
to three dimensions. We always assume a sufficiently large persistence length, such that 
the weakly bending approximation applies and the huctuations perpendicular to the 
mean orientation of the polymer segment are small and can be treated on a Gaussian 
level. 

Our elemental structure is a grafted hlament of contour length L. The probability 
to hnd its tip at position (x, y) with orientation 0, given that it is grafted at (xq, yfi) with 
orientation oj is denoted by ^^.(x, |/, 6'|(xo, |/o, cu). For the simple case of perpendicular 
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Figure 1. Left: Single branched filament grafted on a soft graft and confined in 
the x-direction by a fiat membrane. Thermal bending fluctuations are not shown in 
this drawing. Li is the contour length between graft point and branch point, L 3 is 
the contour length between branch point and filament tip. The contour length of the 
mother filament is L = Li + L 3 , Lb is the branch contour length. The branch angle is 
7 = 70° throughout the study, and w is the tilt angle. The numbers 1, 2 and 3 refer 
to the filament tip coordinates in Eq. 7. Right: Network of branched filaments with 
various orientations and locations of branch points. The length i5o is the distance of 
the filament tips from the graft plane without bending and fluctuations. The lengths 
L and Lb obey L = , Lb = ^ denotes the distance between the graft 

point and the leading edge membrane. 


grafting, a; = 0, at (xo,t/o) = (0,0), Gl obeys [25] 

0/ 

Gl{x, y, 0|O, 0, 0) oc exp[-j^{y^ - LyO + —^)](5(a; - L). (1) 

The general case is obtained by a translation and rotation according to 

9 ^ 9 — u 

y ^ {y — yfi) costu — [x — xq) sintu 

X ^ {x — Xq) coso; + {y — yo) sino; ( 2 ) 

and explicitly given by Eq. 16. 

In a first step we compute the probability, Pt{x\xo), that the endpoint of the tilted 
polymer is at a distance oi x — Xq from the graft point: 

Pt{x\xo) = j j dyd9 GL{x,y,9\xo,t],u}). (3) 

We consider an elastic, fluctuating structure into which the hlament is grafted and hence 

model it by a fluctuating spring in x-direction with spring constant Kg = {kBT)K, 

zero equilibrium length and a distribution of the spring extension proportional to 
exp{—{Kxl)/2). The stiffness of the substrate, K, is assumed to be large as compared 
to the stiffness for transverse fluctuations, K ^ K± = of a hlament of contour 
length L. In this study we take K = 100 3 ^^^ cos{u)p • probability Pfix) to hnd the 
hlament tip at x then follows by averaging over all graft point positions xq-. 

Pt{x) = j Pt{x\xo)exp{-^^). 


( 4 ) 
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The filament exerts a force on an impenetrable flat membrane at a distance 6 from 
the graft plane (see Fig. 1). The force originates from the reduction in the number of 
hlament conhgurations due to the constraint x — xo < S imposed by the membrane; we 
therefore call it an entropic force. The fraction of conhgurations satisfying the constraint 
is given by 

Zt{5) = j dxPt{x )&{6 - x). (5) 

The derivative of Zt with respect to 6 is the entropic force [26] 

fi{L, 6 ,u) = kBT-^lnZfid). (6) 

These results for a single grafted polymer are easily generalised to the branched 
structure shown in Fig. 1. The probability to hnd the two endpoints at X 2 and x^ 
respectively given that the structure is grafted at point xq is calculated by the following 
expression; 



Pb{x2,X3 I Xq) = 

dxidyid9iGLfixi,yi,9i \ xo,0,uj) 
dy2d92GLfix2,y2,92 \ xi,yi,9i - 7 ) 
dy3d93GLfix3,y3,93 \ 


X 



X 



( 7 ) 


The stiffness of the graft is again taken into account by averaging over the positions of 
the graft: 


Pb{x2,X3)(x / dxoPb{x2,X3\xo)exp{- 


Kxl 


) 


( 8 ) 


The explicit expression for Pf,{x 2 ,X 3 ) is given by Eq. 22 in the Appendix. The fraction 
of conhgurations satisfying the constraint obeys in analogy to Eq. 5 


Zb{S) = J dx 2 j dx 3 Pb{x 2 ,X 3 )Q {6 - X 2 )Q{S - X 3 ), (9) 

and the entropic force follows from 

fb{L, Lb, L 3 , 6 , u) = ksT ^ In Zb{L, Lb, L 3 , 6 , uj). (10) 

The force fi depends on the parameters L = (Li + L 3 ), L 3 , oj and S {Lb is hxed by u, 
L and L 3 ). We will present our results as dimensionless quantities, and scale lengths 
by (5o and force by ksT /for that purpose. We comment on which specihc values 
of parameters are suggested by experimental observations in the Discussion and the 
Appendix below. 
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Figure 2. Left panel: force versus relative distance, S/Sq, between membrane and 
graft point for the symmetric case uj = 35° and L 3 = Lb = ^. Right panel: force 
versus tilt angle uj for 5 = Jq. The curves represent the ratio of ^ = 10, 5 and ^ 
from top to bottom in both figures. 


2.2. Properties of single branched filaments 

The entropic force exerted by the branched hlament on the membrane is shown in the 
left panel of Fig. 2 for the symmetric case uj ~ 35° and branching at the midpoint of 
the mother hlament. The force decreases with increasing distance 5 between grafting 
plane and membrane. An inhnitely stiff hlament {Ip very large) would just touch the 
membrane, if 5 equals (5o = L cos uj. But the force exerted by a semihexible hlament 
is nonzero even for <5 > (5o, because the tilted branched structure exhibits huctuations 
with the endpoint reaching beyond 5o- In this regime the force is nearly independent of 
persistence length, whereas for <5 < 5o we observe a strong increase with Ip. 

The force exerted by the hlament on the membrane is crucially ahected by the tilt 
angle. It is plotted as a function of tilt angle uj for hxed <5 = (5o in the right panel of 
Fig. 2. Remarkably, there is a shoulder in all three curves around uj ~ 35°, indicating 
that the symmetric structure generates comparatively large forces. For a; —)■ 0, the 
force would diverge for a stih graft, since Eq. 1 excludes longitudinal huctuations of 
the hlament tip. The high but hnite stihness of the graft limits the force which can 
be exerted on the membrane. These two opposite ehects generate the maximum in the 
force close to w = 0. 

Figure 3 shows the entropic force of a branched structure for the whole range of 
branch point positions on the mother hlament. At the time of Arp2/3 binding to the 
mother hlament (Lb=0), the complete structure has of course the properties of the 
mother hlament. When the branch point has reached the graft plane, the branched 
structure corresponds to two hlaments with the corresponding tilt angles. In between, 
the force has a maximum at Lb = 0.484L. The maximum force is about 2.25 times the 
force of two hlaments with length L. For the special case under consideration, 6 = 60 , 
the maximum of the force scales like the square root of Ip, as shown in the right panel 
of Fig. 3. For the compressed case, 6 < 6 q, we observe a crossover to linear scaling. To 
understand this behaviour of the force, we consider the case of a single hlament: The 
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Figure 3. Left: force versus relative length of the daughter branch. Below: the 
configurations of the branched filament corresponding to LbjL = 0,0.5,1.0. The 
parameter values are: L 3 = Lf,, oj = 35°, 5 = i5o. From high to low force values 
the symbols represent the ratio of ^ = 10, 5 and respectively. Right: log-log plot 
of the maximum force as a function of persistence length /p, fitted to a square root 
dependence f^ax oc ^JTp. 

force has a thermal (entropic) and a nonthermal (energetic) contribution. The thermal 
force dominates for 5 ^ (5o- Evaluating the force given in Eq 21 for 5 = 5o and large K, 
we explicitly see / oc As we compress the filament the energetic force becomes the 
dominant contribution. It can be computed as the force to bend a grafted cantilever 
beam resulting in / = ( • Hence we observe a crossover from the scaling 

of the force with to linear scaling as d is decreased below do 

Both, the maximum force, fmax, as well as its branch point position, , depend 
on the tilt angle u. The general dependence on u is shown in the top left panel 
of Fig. 4 for fmax and in the top right panel for Since the dependencies are 

nonmonotonic, we plot in the bottom panel the configurations which give rise to the 
maximum force and help to understand the non-monotonic dependence. For very small 
u, i.e. almost perpendicular grafting, is very small and hence also the distance 

between the branch point and the endpoint of the mother filament, denoted by 
As the tilt angle, uj, increases, the branch point moves further away from the leading 
edge membrane and hence both, and increase. As u approaches 70°, the 

branch is almost perpendicular to the membrane and the maximum force is obtained 
for branching at the grafting plane, implying = L. Plotting /L as in the top 

right panel of Fig. 4, one actually observes a decrease of because L grows faster 

than as a; —)■ 70°. 

imax jg independent of the persistence length, but does depend on d. For 
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Figure 4. Upper left: Maximum of the force fmax obtained by varying if,, versus 
tilt angle w; we have S = Sq for the curves. From high to low force values the curves 
represent the ratio of ^ = 10, 5 and respectively. Upper right: Relative branch 
length of the maximum force versus tilt angle w for the same parameters 

(independent of the persistence length fi). Bottom: Relative partial contour length of 
the mother filament in between the endpoint and the branch point with maximal force 
j^max _ j^max cos w; also shown are three representative configurations. 


5 > 5q entropic contributions dominate and the force is largest for two independently 
fluctuating filaments, such that /L quickly approaches 1 as 5 extends beyond ho- 

On the other hand, for h < ho) elastic contributions are important. In the symmetric 
case, we hud 0.435 < < 0.484 for 0.5 < h/ho < 1-0. 

The maximum force is observed for ta ~ 0 and cj ~ 70°, because either the mother 
hlament or the daughter hlament are perpendicular to the membrane (upper left panel of 
) . When the branch is perpendicular, the force is mainly determined by the fluctuations 
of the branch point, and when the mother hlament is perpendicular by the huctuations 
of the graft point. Moving away from perpendicular incidence the force has to decrease, 
giving rise to a minimum for intermediate cu. 
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Figure 5. Left panel: Average force per branched structure of the network as a 
function of S/Sq and Ip/Sq. Right panel: The force constant fcj, which is the derivative 
of fb^o/ksT with respect to 5/5q. 
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Figure 6. Ratio R of the average force per filament of branched to unbranched 
networks versus relative distance 5/5q between the membrane and the graft point; 
^ = 10. The squares refer to the random network with a uniform orientation 
distribution in the range —70° < w < 70° as measured in ref. [27], the circles show 
results with a narrow uniform distribution of orientations in 30° < w < 40°, the dots 
refer to the orientation distribution given by Eq. 13 as measured in ref. [28]. 





2.3. Extension to the F-actin network: Properties of an ensemble of branched filaments 

We would like to obtain an estimate of how branching affects the network properties. 
With the theoretical means set up above and in the spirit of the study, we calculate the 
force as an average across an ensemble of branched structures, in which all interactions 
between the branched hlaments are neglected. The ensemble is described by branch 
point (L 3 )and orientation (cu) distributions. The branch point is equally likely anywhere 
on the mother hlament, corresponding to 0 < L 3 < (5o/ coscu. The length of the branch 
obeys Lh = As far as the orientation of the mother hlaments is concerned, 

several scenarios have been discussed in the literature. In ref. [27], electron microscopy 
was used to determine the orientation of hlaments in lamellipodial actin networks. The 
distribution was found to be approximately uniform in the range of angles between 0 
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and 60° with a small contribution between between 60° and 75°. We describe it here as 
a uniform distribution between -70° and 70°. 

The force for hxed (5o averaged over branch point positions and tilt angles is given 
by 


< m >= 


doj 


, ^0 
cos ui 


X 


dU h{- 


coscn 

Ls cos oj 


( 11 ) 


■, L 3 , (5,0;). 


Toso;’ cos(o; — 7 ) 

where Ao; = Umax — (^min denotes the range of the distribution. The average force is a 
monotonic function of both, the persistence length and the distance 6 between grafting 
plane and membrane, as shown in Fig. 5. 

We can compare our results with measurements of the force-velocity relation of hsh 
keratocyte lamellipodia. The forces exerted by the leading edge of the freely running 
cell immediately upon collision with the cantilever were below the force resolution of 
the cantilevers [13, 11, 12], but caused an immediate decrease of leading edge velocity 
by 1-3 orders of magnitude. Hence, the leading edge is much softer than the softest 
cantilever used in the experiments, which had a force constant of 9.1 nN/pm. If the 
leading edge had the same force constant, a single branched structure would need to 
have a constant of about 11.8 pN/pm with 220 hlaments/pm [29] and a contact length 
of about 7 pm [12]. The value of 6q in the freely running cell was estimated to be 
~1.3 pm [13], and such a value is also supported by the data in ref. [30]. Hence, if 
the leading edge had the same force constant as the softest cantilever, the dimensionless 
force constant of a single branched structure would be 11.8 pNpm“^(5o(kBT)“^=4748 
(kBT=4.2 10“^ pNpm). The values in the right panel of Fig. 5 show that the branched 
hlaments are much softer on average, which is in agreement with the dramatic velocity 
drop of the lamellipodium leading edge upon collision with the cantilever. 

We do not know measurements of the pressure exerted by the hlaments on the 
leading edge membrane in the freely running hsh keratocyte but can estimate it from 
our results. Assuming (5o=1.3 pm, Ip/So ~ 10, (5/(5o ~ 0.8 [13] we obtain fdo/ksT ^ 20 
and a force per branched structure of 0.063 pN. The pressure exerted by the hlaments 
on the leading edge membrane is in the range of 46 Pa (with 110 branched structures 
per micrometer and a lamellipodium height of 150 nm as in [12]). 

The value of (5o decreases during the force-velocity measurement due to a dynamic 
equilibrium between polymerization and cross-linking [13]. In the stalled state, 
(5o ~0.27 pm applies (Fig. S3 of [13]). The measured stall pressure exerted on the 
leading edge by the cantilevers is 300-750 Pa [13] and 110-430 Pa [12]. These pressures 
correspond to stall forces of 0.045-0.1125 nN [13] and 0.0165-0.0645 nN [12] per 
micrometer leading edge. The values of fdo/ksT resulting from these force densities 
are in the range 10-70. That entails 6/6q < 1 (see Fig. 5) in agreement with modelling 
results in ref. [13]. 

To assess the effects of branching on network properties, we calculate the ratio R 
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of the force exerted by a branched network to the force of an unbranched network with 
the same angular distribution. We use the same number of hlaments in the unbranched 
network as there are branched structures in the branched network: 

<fb> 


m = 




C^max 1 X / r / r \ ' 

J , duj ft{do cos Uj,d,uj) 

•J U/rn.i.rJ. 


( 12 ) 


The ensemble averages of branched and unbranched networks with the uniform 
orientation distribution between -70° and 70° behave very similar under compression 
(Fig. 6, squares). This ratio has the remarkable property of depending only very weakly 
on 6 . It is almost independent of the persistence length Ip and the graft stiffness K 
as well (data not shown). Hence, the behavior of both networks scales very similar 
in dependence on these parameters. The ratio is about 2, i.e., the average force per 
hlament tip for a given value of 6 is the same for branched and unbranched networks. 

However, there are obvious qualitative differences between single hlaments and 
branched structures illustrated by the non-trivial dependency of the properties of 
branched structures on the branch point position in Figs. 3 and 4. Indeed, if we 
use a narrow orientation distribution around u = 35° the value of R increases 
with increasing compression. Therefore, we also investigate non-uniform measured 
lamellipodial orientation distributions to investigate whether the dependency of R on 
the distribution affects the behavior of lamellipodial networks. Distributions peaked 
either at ca = 0 or ca = ±35° have been observed in refs. [28, 31, 32]. The distribution 
in ref. [28], their Fig. 4, is typical for the measured non-uniform distributions and can 
be approximated by 


P{u) = 0.008012 


e 2-20.5^ -|- e 2-20.5^ 


(13) 


+ 0,0006907, -85° < w < 85° 


The dots in Fig. 6 show the results for R. They are very similar to the results for 
the uniform distribution. Obviously, the width of 40°-45° of the peaks in lamellipodial 
orientation distributions is too large for an essential effect of branching on the parameter 
dependencies. 


3. Discussion 

We investigated the properties of branched hlaments grafted into an elastic graft. Their 
stiffness has a maximum in its dependence on the branch point position. Branched 
structures with the optimal graft point position can be more than four times as stiff as 
a single hlament with the same tilt angle as the mother hlament (Fig. 3), while requiring 
only 1.5 times the polymer length of the mother hlament. 

The mechanical properties of branched F-actin networks depend on their orientation 
distribution. With measured lamellipodial distributions, networks of branched 
structures are about twice as stih as unbranched networks. The total network force 
of branched networks scales essentially the same as the one of unbranched networks 
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with the parameters F-actin persistence length, graft stiffness and compression (6). 
An intnitive explanation wonld be, that in the end it is the single hlament behavionr 
determining the stiffness for both nnbranched and branched hlaments, since the branch 
leans on the (single) mother filament when experiencing a force. 

Onr theory considers individnal branched filaments and thns establishes the 
constitntive relations on which complex network stndies inclnding cross-linking can 
be based. We considered F-actin networks as dehned by their geometrical property 
distribntions withont interactions of hlaments by cross-linking or entanglement. This 
implies that we cannot acconnt for visco-elastic properties. A variety of evidences 
snggests the existence of a region close to the leading edge, where cross-linking is 
not dominating the network properties and to which onr theory directly applies. 
Measnrements of the ratio of nnmber of the cross-linkers to the nnmber of actin molecules 
in hbroblasts show the existence of a gradient for a-actinin and ABP-280/hlamin. 
The number ratios are low in a region juxtaposed to the leading edge with a depth 
of about 1.5 (xm (see Fig. 5 of ref. [15]). Svitkina and Borisy conclude from these 
results and structural information from electron micrographs that the impact of ABP- 
280 and a-actinin on hlament cross-linking is likely to be expressed more deeply in the 
cytoplasm [15]. Their statement is supported by the immediate response of the leading 
edge of hsh keratocytes to small forces indicating also weak cross-linking in the network 
region close to the leading edge [13, 11, 12]. This suggests that while understanding 
of the visco-elastic properties of the network in the lamellipodium bulk requires taking 
cross-linking into account, our ensemble average is applicable to a network region close 
to the leading edge. The reproduction of both the weak and strong force responses 
of the lamellipodium leading edge measured in force-velocity relations by our network 
calculations strongly supports that conclusion (Fig. 5). 

We did not take contributions from entanglement or excluded volume ehects into 
account when calculating the network forces. This implies that R provides only a 
meaningful approximation, if these ehects are similar in branched and nnbranched 
networks. To the best of our knowledge, that has not been investigated quantitatively 
yet. We can only provide heuristic considerations in favour of our assumption based 
on comparing a variety of simulations with and without excluded volume ehects with 
force-velocity measurements. 

Model networks of semi-hexible hlaments not taking into account excluded 
volume ehects reproduce the elastic properties measured in force-velocity experiments 
quantitatively [13]. Schreiber et ah simulated the force-velocity relation of motile cells 
with rigid rods as model hlaments taking excluded volume ehects into account [33]. They 
found excluded volume ehects to be stronger in branched than in nnbranched systems. 
The model network of Schreiber et ah shows a response to external forces in the force- 
velocity relation at about 8 nN/qm [33]. However, the lamellipodium leading edge 
exhibits elastic responses to forces smaller than 0.05 nN/qm in experiments [13, 11, 12]. 
Additionally, bending of hlaments has been observed in lamellipodia [27, 34, 35], i.e., 
hlaments do not behave like stih rods. Hence, the lamellipodium network is likely to be 
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in a parameter regime where excluded volume effects are less relevant than suggested 
by a network of stiff rods. 

Branching has also been observed with microtubule [36]. The branching angle 
varies between 0 ° and 90 and it is not known, how rigid the connection to the mother 
hlament is. If the branch is rigidly connected, our results should apply also to these 
branched structures with the adapted persistence length (a few millimeters [37]) and 
7 -values. 

In summary, single branched hlaments with intermediate branch point positions and 
tilt angles exert larger forces and are stiffer than two unbranched hlaments. The stiffness 
of a whole network of branched hlaments is largest, if the orientation distribution is 
sharply peaked around ± 7/2 = ±35°. For lamellipodial orientation distributions, the 
stihness of branched and unbranched networks scales approximately the same with a 
variety of parameters (Fig. 6 ), suggesting that the ehect of branching on network stihness 
can be accounted for by rescaling the hlament number of an unbranched network. These 
results are in agreement with the elastic properties of lamellipodia found in force-velocity 
measurements with hsh keratocytes [13, 11, 12]. 


4. APPENDIX 


4 . 1 . Tilted filament 

The Green function of a free semihexible polymer in two spatial dimensions (see Fig. 1) 
and in the weakly bending limit satishes the partial diherential equation (PDF) 


1 ^ + ^^ 
5s dy 


1 52 

TpW 


]G{s,y,e I 0, 1 / 0 , 0 ;) = 0. 


(14) 


The arc length of the hlament contour is denoted s here. The boundary condition 

\iinG{s,ys,0s\O,yo,uj) = S{e-uj)6{ys - yo) (15) 

s—>-0 

realizes the graft. The solution is Eq. 1, after a switch from the coordinates {s,y,6) to 
{x, y, 6). We obtain the Green function for a tilted hlament simply by a rotation: 

GLix,y,e\xo,yo,u) (X (16) 

31 

exp[-j^{{y - yo) cos{u) - (x - xq) sin(o;))^ 


Ir 


31 

X exp[+j^{{y - yo) cos{u) - (x - xo) sin(u;))(0 - o;)] 

X (5[(a: — Xo) cos(o;) + (y — yo) sin(o;) — L]. 

The probability distribution of the position x of the endpoint follows by integration (see 
Eq. 4) 

{x — Xo — L cos(a;))^. 


(J: 


to 


Pt{x\xo) oc exp( 


(17) 
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where 


4L3 


sin cu 




3L 


(18) 


Representing the stiffness of the graft by a spring and averaging over all grafting points, 
Xq, yields: 




Nt - ' 

Here acconnts for the proper normalization and 

2 

^2 _ ^2 I ^ 

The entropic force on the wall is then explicitly given by: 

/ (5—Lcos(lj))‘^\ 

2A:bT exp(-^- 

Jt{d) = 


(\/T(Tf) erfc( 


L cos(lj)—5 > 


(19) 

( 20 ) 

( 21 ) 


4 . 2 . Branched filament 


Starting from Eq. 7 and averaging over graft point positions xq, we obtain 


T-, / ^ /(det(M ))2 

Pb{x2,X3) = (-- )exp{r]iMijr]j), 


TT 


where 


T] = 


^X 2 — Li cos (a;) — Lfe cos(ci; — 7 ) 
X 3 - (Li + L 3 ) cos(a;) 
and M is a 2 X 2-Matrix with components 


Mil — + -^3)^sin(a;)^ -|- -Ip] 

Mi 2 = M 21 = -^[KL\{Li -\- -L 3 ) sin(a;)^ -|- -Ip + 

3 

-KLiLb{Li 2 L 3 ) sin(a;) sin(a; - 7 )] 


M22 — -^[KL‘ 1 {?)Li + Lb) sin(ci; — 7)^ -|- -Ip + 
3 KLlLbSm{uj) sin(a; — 7 ) + KLl sin(a;)^], 

C = R'sin(a;)^[(Li -|- ^L3)LfL3 sin(a;)^ 

o 

— 2(Li — 2Ll)L\LbL3 sin(ci;) sin(ci; — 7)] 

-h LIK sin(a;)^ sin((:n - ^)‘^[L‘\ + ^{L\Lb + LbL\) 

o 

+ ALlLbL 3 + AL^Ll{Lb + L 3 )] 

-|- 6/p[(Li -|- -L3)LI sin(ci;)^ — 2LiLbL3 sin(ci;) sin(a; — 7) 
o 

+ (^1 + g^6)7 sin(w — 7)^] 


( 22 ) 

(23) 

(24) 
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4 . 3 . Parameter values 

The independent parameters of the model are the F-actin persistence length Ip, the 
equilibrium distance of the hlament tips from the graft plane (5o, the stiffness K of 
the graft, and the branching angle 7 which is hxed by the Arp2/3 complex at about 
70° [38, 39, 40]. The branch point position is uniformly distributed and the tilt angle 
u according to the above distributions. The tilt angle and (5o fix the mother hlament 
contour length as L = 60 / cos(a;). 

We explained in the Discussion, that a region juxtaposed to the leading edge 
membrane with a width of 1.0-1.5 pm is similar to an experimental realization of the 
network conhguration in Fig. 1, since it is weakly cross-linked. Values for (5o suggested by 
these observations are in the range 1.0-1.5 pm. Branching will have an effect on network 
properties in a conhguration like shown in Fig. 1, if most hlaments are branched, i.e., if 
do is larger than the average branch distance. 

The average branch distance in steadily moving cells has been a matter of debate 
in recent years. Svitkina et ah concluded 20-50 nm from early electron micrographs 
of lamellipodia from Xenopus keratocytes and hbroblasts [15], and 50-200 nm from 
another study in hbroblasts [41]. Later measurements substantially increased that 
value. The average branch distance has been determined to be about 800 nm for 
B16 melanoma cells and hsh keratocytes in ref. [35]. Other studies provide number 
densities of branch points per lamellipodium area for 3T3 hbroblasts. Calculating the 
average branch distance from that branch point density implies assumptions on the F- 
actin concentration. The number of hlaments per micrometer lamellipodium width in 
a distance of 0.1-1 pm from the leading edge is 170-190 [29]. Taking into account that 
hlament orientation is approximately uniformly distributed between 0° and 60° [27], 
this density means 1.16(170-190) pm hlament contour length per pm^ lamellipodium 
area. The factor 1.16 arises from averaging over all tilt angles of the hlaments. Yang 
and Svitkina measured 277 branch points/pm^ [40] in the same sample, i.e. an average 
branch point distance in terms of contour length between 700 nm and 800 nm. Small 
et ah measured less than 225 branch points/pm^ [30], i.e. an average distance of more 
than 860 nm. 

The in vivo persistence length is not known. We can only conclude a reasonable 
range from in vitro measurements. Results from huctuation analysis and measurements 
of network elastic properties yield values of the in vitro persistence length of 15- 
18 pm [37, 42, 9, 43]. The hlaments in these experiments were stabilized with phaloidin 
which most likely increases the value of Ip. Filaments labeled with rhodamine but not 
stabilized with phaloidin exhibited values of Ip between 9 pm and 13.5 pm [43]. The in 
vivo persistence length might be even shorter, since cohlin can substantially reduce it 
even down to 2.2 pm [44, 45]. 
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